epl draft 



Doping evolution of itinerant magnetic excitations in Fe-based 
^oxypnictides 

o: 
o- 

(NM.M. Korshunov 1 ' 2 and I. Eremin 1 - 3 

~^> . Max-Planck-Institut fur Physik Komplexer Systeme - 01187 Dresden, Germany 

L.V. Kirensky Institute of Physics, Siberian Branch of Russian Academy of Sciences - 660036 Krasnoyarsk, Russia 
^3 . Institute fur Mathematische und Theoretische Physik, TU- Braunschweig - 38106 Braunschweig, Germany 

(N . 

ll 
O 
o 

u 

Oh 



c 
o 
o 



> 

en 

o 

00 

O 



PACS 74 . 20 . -z - Theories and models of superconducting state 
PACS 74 . 25 . Ha - Magnetic properties 
PACS 75 . 30 . Fv - Spin-density waves 

Abstract. - Employing the four-band tight-binding model we study theoretically the doping 
dependence of the spin response in the normal state of novel Fe-based pnictide superconductors. 
We show that the commensurate spin density wave (SDW) transition that arises due to interband 
scattering between the hole a-pockets and the electron /3-pockets disappears already at the doping 
concentration x ~ 0.04 reflecting the evolution of the Fermi surfaces. Correspondingly, with 
further increase of the doping the antiferromagnctic fluctuations are suppressed for x > 0.1 and 
the Imx(QAFM,w) becomes nearly temperature independent. At the same time, we observe that 
the uniform susceptibility deviates from the Pauli-like behavior and is increasing with increasing 
temperature reflecting the activation processes for the a-Fermi surfaces up to temperatures of 
about T = 800K. With increase of the doping the absolute value of the uniform susceptibility 
lowers and its temperature dependence changes. In particular, it is a constant at low temperatures 
and then decreases with increasing temperature. We discuss our results in a context of recent 
experimental data. 
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Introduction. — The recent discovery of super- 
conductivity in the iron-based layered superconductor 
La(0i_2;F x FeAs) with T c w 26K [1] has generated the re- 
newed interest in high-temperature superconductivity due 
to consequent development of materials with higher T c 's 
up to 55K that contain other rare-earth elements such as 
Ce, Nd, Sm [2-4] instead of La. The physical properties 
are considered to be highly two-dimensional; the crystal 
structure is tetragonal and consists of the LaO and the 
FeAs layers which are stacked along the c-axis. Similar 
to many layered transition metal oxides the superconduc- 
tivity in oxypnictides occurs upon introducing doping of 
cither electrons [1-4] or holes [5] into the FeAs layers and 
the parent material shows antiferromagnctic transition at 
around 150K [1,6-9]. At the same time, in contrast to lay- 
ered cuprates the parent material remains a metal. The 
observed magnetic moment per Fe atom has been reported 
to range between 0.25/xb [9] and 0.36/j.s [7] and lies in the 
afe-plane. 

There have been various proposals to explain the ori- 
gin of antifcrromagnetism in these systems. Recent theo- 



retical studies suggest several different explanations vary- 
ing from LaOFcAs being an antiferromagnctic scmimctal 
[11-13], or the system with frustrated magnetic ground 
state with two interpenetrating antiferromagnetic square 
sublattice [14-17]. However, the resulting magnetic mo- 
ments is larger than that found in experiment thus requir- 
ing an inclusion of strong fluctuations effects that would 
reduce the magnetic moment. At the same time, starting 
from the purely itinerant models it has been also proposed 
that LaOFcAs has an antiferromagnctic spin density wave 
instability due to the interband nesting of the electron and 
the hole Fermi surfaces [6,9,18]. The resulting magnetic 
moment has been found to be about 0.33/is which agrees 
with experimental data. Despite the right order of mag- 
nitude it remains to be seen whether the strong electronic 
correlations that might be important in LaOFeAs due to 
the Hund exchange [19,20] will modify this result. It has 
been also argued that a combined effect of spin-orbit cou- 
pling, monoclinic distortions, and p — d hybridization may 
invalidate the simple Hund coupling scheme [21]. 

In order to understand how the magnetism and the spin 
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fluctuations in La(Oi_ K Fa;FeAs) evolve as a function of 
doping in this letter we present the study of the mag- 
netic excitations using the tight-binding scheme adopted 
previously [22]. In particular, we show that the commen- 
surate spin density wave (SDW) transition that arises due 
to interband scattering between the hole a-pockets and 
the electron /3-pockcts at the Fermi surface (FS) disap- 
pears already at the doping concentration x at 0.04. Cor- 
respondingly, with further increase of the doping the anti- 
ferromagnetic fluctuations arc suppressed and at x k, 0.1 
the Ihix(Qafm,w)/w becomes nearly temperature inde- 
pendent. At the same time, we observe the uniform sus- 
ceptibility deviates from the Pauli-likc behavior and is in- 
creasing with increasing temperature reflecting the activa- 
tion processes for the a-Fermi surfaces up to temperatures 
of about 800K. With increase of the doping the absolute 
value of the uniform susceptibility is decreasing and its 
temperature dependence changes. It is a constant at low 
temperatures and then decreases with increasing temper- 
ature. 

Theory. — The effective low-energy band structure of 
the undoped LaOFcAs can be modeled by the following 
single-electron model Hamiltonian for the folded Brillouin 
Zone (BZ) with two i^e-ions per unit cell [22]: 



Ho = - 



k,Q,(T 



e n^ta 



k,?'.(7 



k"k?fj^ k ^' 



(i) 



where i = 
€ l are the 

j.ai,a2 



011,012, 01, 02 refer to the band indices, 
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electron /3-pockcts around the M point of the folded 
BZ. Using the abbreviation (e l ,t\,t 2 ) we choose the pa- 
rameters (-0.60,0.30,0.24) and (-0.40,0.20,0.24) for the 
cm and «2 bands, respectively, and (1.70,1.14,0.74) and 
(1.70,1.14,-0.64) for the 0\ and 02 bands, correspond- 
ingly (all values are in eV). 

In Fig. 1(a) we show the resulting energy dispersion 
along the main symmetry directions of the first BZ for 
the undoped case, x = 0. The band structure parame- 
ters were chosen to correctly reproduce the LDA Fermi 
surface topology and the values of the Fermi velocities for 
the hole a and the electron pockets. In particular, we 
have selected the on-site energies and the hopping matrix 
elements assuming the compensated metal at x = and 
calculating the chemical potential self-consistently for the 
filling factor n — 4 (we further assume that there exists 
another band below the Fermi level which is fully occu- 
pied and not considered here). As a consequence, the hole 
Fermi surfaces shifted by vector (jr, n) is nearly completely 
nested with that of the electron pockets in full agreement 
with ab initio density functional calculations [6,19,23,24]. 
Additionally, we take into account the details of the elec- 
tronic dispersions of the bands forming the Fermi surface 
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Fig. 1: (a) Calculated energy dispersion along the main sym- 
metry points of the first BZ for the undoped, x = 0, case, 
(b) LDA band structure (green crosses) and ten-band model 
dispersion (red curves) after K. Kuroki et al. [24]. Note the 
difference in BZ directions in (a) and (b). The large arrows 
indicate the points where bands cross the Fermi level. 



pockets. In order to visualize the comparison, in Fig. 1(b) 
we present the LDA band structure and the realistic ten- 
band model dispersion from Ref. [24]. Note, for other dop- 
ing concentrations the position of the chemical potential 
was deduced from the equation n = 4 + x. 

The resulting doping-dependence of the physical sus- 
ceptibility as obtained by the sum of all interband and 
intraband susceptibilities is shown in Fig. 2. For the un- 
doped case our results are in qualitative agreement with 
that of K. Kuroki et al. [24] and S. Raghu et al. [25], and 
with the ab initio results of J. Dong et al. [6]. 

Within random phase approximation (RPA) the spin 
response has a matrix form: 



XRPA(<1, '^r, 



rxo(q,i^™)] 1 xo(q,iw,„) (2) 



where I is a unit matrix and xo(q, i^m) is 4 x 4 matrix 
formed by the interband and the intraband bare suscep- 
tibilities. For the four-band model considered here the 
effective interaction consist of the on-site Hubbard intra- 
band repulsion U and the Hund's coupling J. There is 
also an interband Hubbard repulsion U' , which, however, 
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Fig. 2: Calculated doping dependence of the real part of the 



physical spin susceptibility, Xo(q, 0) 
refer to the band index. 



Ei;Xo J (q-°). where 
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does not contribute to the RPA susceptibility. The vertex 
is given by 



U J/2 J/2 J/2 

J/2 U J/2 J/2 

J/2 J/2 U J/2 

J/2 J/2 J/2 U 
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For the given Fermi surface topology the main magnetic 
instability in the folded BZ occurs at the antifcrromag- 
netic wave vector Qafm = (t, tt) due to the interband 
nesting between the hole a- and the electron /3-bands 
[6,22,23,25]. This is also clearly visible from our Fig. 2. 
Note that in the unfolded BZ with one Fe-ion per unit 
cell, the wave vector is Q! afm = ( n 7 0) which corresponds 
to the 'stripe'-like ordering of the Fe-spins as observed 
by neutron scattering [7]. Setting the Hund's coupling 
to J = 70meV and choosing U = 320mcV we obtain 
the ordering temperature T/v = 138K as determined by 
det [I — rxo(q,iui m )] — 0. Solving the condition for the 
SDW instability below Tjy which can be regarded as a 
mean-field equation for the SDW order parameter, we ob- 
tain Asdw(T = OK) = 31meV which corresponds to the 
magnetic moment per two Fe sites to be \i « 0.33/xs. 

Note, the small values of U and J used here is a conse- 
quence of the absence of the self-energy corrections within 
RPA approach. Such corrections would reduce the value of 
the absolute magnitude of the spin susceptibility and cor- 
respondingly yield larger values of the coupling constants 
U and J. 

In the inset of Fig. 3(a) we show the doping dependence 
of the Neel temperature. One finds that it decreases quite 
rapidly as a function of doping and goes to zero already 
at x ra 0.04. The reason of the rapid suppression of the 
Neel temperature is quite obvious. Away from x = the 
spectral weight of the hole a-pockets at the Fermi sur- 
face is decreasing and the condition for interband nesting 
becomes worse as it is readily seen from Fig. 2. In particu- 
lar, one finds that the peak at the antifcrromagnctic wave 
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Fig. 3: Calculated lim^-,0 ImxRPA(QAFM,uj)/u for various 
doping concentrations. Note the log scale in (a). The inset in 
Fig. 2(a) shows the calculated doping dependence of the Neel 
transition temperature. 



vector, Qafm, decreases quite rapidly away from x = 0. 
Remarkable that such a small deviation from undoped case 
changes the situation dramatically also in the NMR exper- 
imental data [10] that gives an additional support in favor 
of the nesting scenario of the antiferromagnetic transition. 

In Fig. 3 we show the doping dependent evolution of 
the lim^^o Imx(Q>iFAf> u>)/w. One finds that at x = 
it diverges at Tjv and with further increase of the doping 
the antiferromagnetic fluctuations are quickly suppressed. 
Remarkable one finds that hnx(QAFM , w)/w at x = 0.12 
docs not show any enhancement characteristic for strong 
antiferromagnetic spin fluctuations and stays nearly con- 
stant. Despite the fact that the real part still shows the 
peaks around Qafm the antiferromagnetic fluctuations 
are quite strongly suppressed in the imaginary part of the 
spin susceptibility. This is due to the fact that the RPA 
response has a matrix form and thus the damping of the 
fluctuations is quite strong. At the same time, one has 
to keep in mind that the spin dynamics on the As sites 
originating from the stripe-like ordering of the Fe spins as 
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Fig. 4: Calculated uniform part of the total spin susceptibility, 
Xrpa(<\ — * 0,w — > 0) as a function of temperature at x=0 
(a) and for other doping concentrations (b). In (a) we show 
also the partial contributions including the interband and the 
intraband transitions. 



probed by NMR will be suppressed due to hyperfine inter- 
action. Therefore, further experimental studies are neces- 
sary to understand the evolution of the antiferromagnetic 
fluctuations in these systems. 

In Fig. 4 we show the temperature dependence of the 
uniform susceptibility for various doping concentrations. 
It is interesting to note that the uniform susceptibility 
above T^ does not show either the Pauli-like or the Curie- 
Weiss type behavior. At zero doping concentration the 
total susceptibility is increasing as a function of tempera- 
ture up to 600K and then decreases following the Curie- 
Weiss like behavior. Looking on the partial contributions, 
one finds that this temperature dependence is determined 
mainly by the transitions between a-bands which produce 
the hole pockets around the r-point in the BZ. In particu- 
lar, the hole-like Fermi surfaces of the a± and a2 bands are 
only slightly splittcd. The gap between the two Fermi en- 
ergies that would occur for zero transferred momentum is 
about 50meV. Therefore, due to the temperature activated 
transitions between two a bands, the interband suscepti- 
bility will increase with increasing temperature up to 600K 



and then decrease. This slight increase of susceptibility is 
quantitatively consistent with available experimental mea- 
surements [9,26]. Note that the transitions within (3 bands 
and between a and (3 bands show almost Pauli-like behav- 
ior. Upon changing doping the overall magnitude of the 
susceptibility is decreasing which reflects the reduction of 
the total susceptibility as also shown in Fig. 2. The lat- 
ter occurs due to the filling of the hole pockets. We also 
observe the change in the temperature dependence of the 
uniform susceptibility. For x > 0.1 the uniform suscep- 
tibility is constant up to 200K and then decreases with 
increasing temperature. This change occurs due to the 
filling of the a-bands upon varying doping and the reduc- 
tion of their relative splitting as can be seen from Fig. 1. 

We finally note that in our analysis we assume all ma- 
trix elements for calculations of the spin susceptibility to 
be unity. Note that a qualitative agreement between our 
results and those found in Refs. [6,24,25] for the undoped 
case seems to justify our approach. In addition we further 
neglect the other three-dimensional band that gets quickly 
filled upon doping. Although its inclusion may be im- 
portant with regard to the formation of three-dimensional 
Neel order, it will not change much the doping dependence 
of the two-dimensional in-plane magnetic fluctuations. 



Conclusion. — We have analyzed the doping depen- 
dence of the spin excitations in La(Oi_ K Fa;)FeAs based 
on a purely itinerant model. We find that the interband 
antiferromagnetic spin fluctuations are rather rapidly sup- 
pressed and the Neel temperature disappears already for 
x w 0.04. With further increase of the doping the short- 
range antiferromagnetic fluctuations disappear at a; w 0.12 
in agreement with NMR data. Given the fact that the 
superconductivity seems to be strongest at this doping 
concentration it is interesting to see whether these fluctu- 
ations can be responsible for the formation of supercon- 
ductivity. Due to the multi-orbital character the uniform 
susceptibility shows neither Pauli-like neither Curie-like 
behavior. In particular, for small doping the total sus- 
ceptibility is increasing up to 600K and then decreases. 
With increasing doping the susceptibility stays constant 
at small temperatures and then lowers with increasing 
temperatures. We find that this characteristic behavior 
originates from the interband transitions between slightly 
splittcd a-bands. 

Note added: After submission of this manuscript we be- 
came aware of the study by Anisimov et al. [27], where the 
values of the average Coulomb repulsion U and Hund's ex- 
change J were obtained by the first principles constrained 
density functional theory in Wannier functions formalism. 
Due to the derealization of Wannier functions for the 
Fe — 3d basis set, the Coulomb parameters were signif- 
icantly reduced in comparison to their intraatomic values 
and became U « 0.6eV and J ps 0.5eV. These are close to 
the effective parameters used in the present study. 
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